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Abstract 

Explicit representations of densities for linear parabolic partial dif- 
ferential equations are useful in order to design computation schemes 
of high accuracy for a considerable class of diffusion models. Approx- 
imations of lower order based on the WKB-expansion have been used 
in order to compute Greeks in standard models of the interest rate 
market (cf. However, it turns out that for higher order approx- 

imations another related expansion leads to more accurate schemes. 
We compute a local explicit formula for a class of parabolic problems 
and determine a lower bound of the time horizon where it holds (given 
a certain bounded domain). Although the local analytic expansions 
hold only for strictly elliptic equations we show that the expansions 
can be used in order to design higher order schemes for various types 
of (micro)-hypoelliptic and semi-elliptic equations, e.g. the reduced 
market models considered in [7] or front fixing schemes for multivari- 
ate American derivatives [S]. 

2000 AMS subject classification: 60H10, 62G07, 65C05 

1 Introduction 

Higher order approximations of fundamental solutions or densities for scalar 
linear partial differential equations of parabolic type are very useful as an 
element for the design of higher order schemes and efficient algorithms in 
the context of a large class of diffusion models. Note that this concerns also 
models which do not satisfy the rather strict conditions which we need in 
order to get local convergent expansions of densities for such equations. As 
examples consider the front-fixing iteration for American options in |3j or 
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the analytic AD-scheme in [7] for semi-elliptic equations of type 



du 

dt ~ 



1 v^m. 4 2 

2 2^i=l 



U + AqU 



(1) 



u{Q,x) = fix) 



on the domain [0,T] x M", and where 



n 



d 




(2) 



and are smooth vector fields vector fields which satisfy satisfy the Hormander 
condition only with respect to the subspace of dimension d < n at each point 
X, i.e. the sets 



{A,[Aj,Ak],[[Aj,Ak],Ai],---\l<i<m, < j,k,l ■ ■ ■ < m} (3) 



span a linear subspace Wx of dimension d < n at each point x G M". If 
d = n then the equation ([1]) is micro-hypoelliptic and therefore hypoelliptic. 
Here [., .] denotes the Lie bracket of vector fields. Furthermore, recall that a 
differential operator L with C°°-coefficients is called hypoelliptic on an open 
set C if Lu G C°° implies u in C°° for any distribution u. Operators of 
the type of equation ([1]) such that ([3]) spans the full space for all x E M"' 
, are microhypoelliptic, i.e. preserve the wave front set 



where WF(it) is the intersection of the characteristic varieties of pseudo- 
differential operators P of order zero which satisfy Pu G C°° . Well, it 
is clear that the operators of ([T]) are not microhypoelliptic for d < n \n 
general. This is reflected by the fact for pseudodifferential operators of 
negative order the characteristic variety equal the whole co-tangential bundle 
of the domain Q.. Hence, the microhypoelliptic theory is designed for 
pseudodifferential operators of order (at least), and this does not apply 
in general if the Hormander condition is satisfied on a subspace of lower 
dimension than n. 

Note that typical diffusion models in finance, e.g. (multivariate versions) 
of the Heston model, are often microhypoelliptic or hypoelliptic, but not 
strictly elliptic. Moreover, in many practical situations of related stochastic 
ordinary differential equations the number of Brownian motions may be 
less than than the dimension of the domain. The situation is considered in 
[7]. Another example of semi-elliptic problems which are not hypoelliptic 
in general, are free boundary problems. Consider for example American 
derivatives in finance. Let WJ^ denote the set of m-tuples of strictly positive 
real numbers. Then for each starting point (i, x, y) G [0, T] x O C x 



WF(u) = WF(Ln) 



(4) 
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a typical diffusion market model is described by the stochastic differential 
equation 

(5f'^y/'^) = (x,y)G0, 

^ = r{s, Sl^^)ds + E •=! <^k,{s, St'\ Ys'^dWi, (5) 

which has a unique strong solution, if the drift functions 
r : [0, T] X ^ M, 

(6) 

zv; : [0, T] X M'^ ^ M, / = I,-- - 
and the (volatility of) volatility functions 

A,- : [o,r] xM'^^M, I ■■,?}, iG{i,---,4 (7) 

satisfy certain Lipschitz conditions. In this context the value function V{t, x, y) 
of an American option with payoff <j) is given by 

V{t,x,y) := sup Eq 
reStopfjy, 

where Stopj^ set of all stopping times with value in [t, T]. The related 

obstacle problem is semi-elliptic. Such a problem is (micro)-hypoelliptic only 
on a subspace in general. Next a natural question occurs: can we use approx- 
imations of fundamental solutions of strictly elliptic equations, i.e. equations 
of form ([T|) which satisfy a uniform ellipticity condition on the whole space, 
in order to design algorithm in situations of (micro)-hypoellipticity or semi- 
ellipticity? More precisely and more general, assuming regularity and at 
most linear growth of the coefficients, how may we use local expansions of 
the fundamental solutions p of form 



p{t,x]0,y) = 




in the context of the semi-elliptic types of equations mentioned above? 

Remark 1.1. Note that representations of the form ([9]) hold in the time- 
homogeneous case. In the time-inhomogeneous case they are of the form 



p{t,x;0,y) = 




and they are computerized in a similar way (see below). 
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Here, the dk are solutions of recursively defined linear partial differential 
equations of first order, and dji is a Riemannian metric defined by the line 
element 

n 

ds'^ = a*'^{x)dx,dxj. (11) 

Remark 1.2. Note that expansions of type (|10p are different from the WKB- 
expansion considered in [2] and in [3]. 

Note that the coefficients a**-' denote the components of the inverse of 
the diffusion matrix {a*j) where the latter coefficients are determined by the 
vector filed coefficients aji in the usual way. The interest in such expansions 
is related to the observation that the first order equations which determine 
the coefficient functions (i^ are easier to solve than the original second order 
equations. Well, it is sometimes not that easy. In order to determine the 
line element (|lip we need to solve nonlinear eikonal equations of the form 

1 " 

Furthermore, we have to construct solutions of ()12p in regular spaces since 
the higher derivatives of dj^ appear as coefficients (and are involved in source 
terms) of the equations for the higher order terms d^ in order to get accurate 
data for these equations. Note that the 'iff '-condition on the 'boundary', i.e. 
the condition d^(x,y) = if and only ii x = y, leads to a term of lowest 
order of form 

n 

J2a*'^{y){xi-yi){x,-yj) (13) 

for d^. Prom that term it looks hopeless to search for analytic expansions 
of type ([9]) in the context of (micro)-hypoelliptic models. In such cases the 
points of elliptic degeneracies of such models are related to singularities of 
the inverse (a*-') which defines the line element. However, this does not 
mean that it is impossible to find good approximations of the density in the 
form ([9]) in many cases. Let us consider an example in finance. Univariate 
or multivariate stochastic volatility diffusion models like (0) are usually of 
form 

dS = fi{S)dt + a{S,Y)dW (14) 

with /i some drift vector, and a some dispersion matrix-valued function. Lets 
assume that elliptic degeneracies of aa^ appear for some set of arguments 
of measure zero (this is true for most of the standard stochastic volatility 
diffusion models such as the Heston model and multivariate versions of it). 
Indeed Malliavin calculus tells us that for models satisfying we typically 
have related weakly invertible covariance matrices. More precisely, if n = d 
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and the Hormander condition holds, consider the associated Stratonovic 
integral of a process starting at x G M" is 



Xt = x+ / AoiX,)ds + y] / AkiX,)odWkis). (15) 
Jo f^^^ Jo 

Then the associated covariance matrix process at of the process satisfies 

a^' G LP, (16) 



where 

at = zr' 



[ ZsAs{Xs)Aj{Xs)Z] 
Jo 

and where Z is a matrix-valued invertible process defined by 



Z-''^, (17) 



" j-t 

Zt = h- ZsDAo{X,)ds-y2 ZsDA{Xs)doW,{s). (18) 
Jo Jo 

Here, 1^ denotes the d-dimensional identity matrix. This indicates that the 
set of degeneracies is rather thin (of Lebesgue- measure zero) . This is indeed 
what we observe in finance usually (cf. the Heston model or the SABR- 
model). This leads us to the following consideration. Let D^,^^ be the set of 
arguments S,Y where the lower ellipticity constant of aa^{S,Y) is less or 
equal to e. Since the leading order term of Q, i.e. 

Vexpf-#V (19) 



goes rapidly to zero for arguments in D^^^ if epsilon is small one may consider 
approximations of the density p which are defined to be zero on the set 
D^^g and equal the expansion for some strictly elliptic operator with ellip- 
ticity constant e > in the complementary domain. Well, one has to control 
the time parameter t to be not too large such that p9|) really dominates the 
higher order terms dk where dj^ and derivatives of dj^ are involved. Then 
using the semigroup property one can set up weak higher order schemes (as 
considered in [2] for example) which are time-discretized according to the 
nature of the degeneracies. The analysis of the time-discretization may be 
complicated sometimes, but it is a possible way of approximation in many 
situations of hypoelliptic operators considered by practitioners. 

Prom the computational point of view there are challenges other than 
degeneracies. Let us consider a few. First, the WKB-expansion 

72 



p{t,x,y) = ^^exp ^-^^^ + ^^c,(x,y)t'= j , (20) 
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considered in [2] can lead to numerical instabilities if higher order terms 
Cfc with k > 2 are considered (see discussion below). In [2] the expansion 
was considered up to the first order term ci. The accuracy and efficiency 
of schemes based on analytic expansions is demonstrated by the fact that 
options in the Libor market with ten years of maturity are computable in 
a typical market situation in one time step. However, for higher volatilities 
higher order terms are desirable, and in that case numerical instabilities 
may appear. In order to have a damping leading order term as in ([9]) we 
compute the equations for instead of the in f2U\i . The resulting first 
order equations for dk are more difficult to solve than those for Ck- In this 
paper we compute recursive solutions for the in terms of the c^, and 
recursively explicit solutions for the coefficients dk in the case where there 
is a global transformation of the second order part of the operator to the 
Laplacian (the 'reducible case'). In that case we also derive a lower bound 
for the radius of convergence. The use of first order solutions for the d^ is not 
restricted to semi-elliptic equations with a thin set of elliptic degeneracies. 
For example high-dimensional models in finance are often reduced in order to 
get computationally feasible models. In [T] and in [7] reduced market models 
are considered which lead to semi-elliptic equations which are not micro- 
hypoelliptic (they are micro- hypoelliptic only on a subspace). In general 
such equations do not have regular densities. Indeed they typically have 
densities only in a distributional sense. For example, the lower dimensional 
Cauchy problem on [0, T] x M^: 

du 1 „2 d^u I , , du 

(21) 

^ ^(O,^) = f{xi)+g{x2). 
has a 'distributional density' of the form 

p{t, X, y) := exp {- ^^^2aH^ ^ ^^^'^ + f^t - y2)- (22) 



'2-Kt 



In [7] analytical AD-schemes were defined based on analytical expansions of 
densities considered in this paper. We shall review these schemes below and 
show how analytic density approximations of this paper can be used. In any 
case, the really hard part is the computation of the Riemannian metric and 
its derivatives. One approach based on regular polynomial interpolation is 
considered in [5]. However, we shall see that a more efficient method exists 
based on analytic schemes which we shall consider in a subsequent paper 
(second part of this work). The linear first order equation for the higher 
order terms are sometimes explicitly solvable, especially in the reducible case 
where a global transformation of the second order terms to the Laplacian 
exists. In general they may be computed by the regular interpolation method 
considered in [B]. Our first observation in this paper is that we can obtain 
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recursive expressions for the from recursive expressions for the Ck of the 
WKB-expansion, i.e. in order to compute a local solution of 

ij=l i 

(with analytic data) in the form 



p(t,x;0,y) = ^^=n exp (^-^^ \ Y,dk{t,x,y)t^j . (24) 



we may first compute the WKB-expansion 

1 / dl 



p{t,x- 0, y) = -= exp - -J + J]] Cfc (t , X, y)t^ . (25) 



k=0 



Indeed the connection between the Ck and the dk can be easily computed 
using Leibniz rule. We have 

do = exp(co), (26) 

and 

k . 

dk = '^^dk^iCi. (27) 

i=l 

The case where the diffusion part can be globally transformed to a Laplacian 
is of special interest, because this makes it easier to study the convergence 
behavior of the higher order terms dk and the horizon of convergence. We 
call this case the reducible case. Accordingly, in the first part of this paper 
we derive explicit analytic formulas of fundamental solutions to scalar linear 
equations of the form 

- = An + ;^Mt,x)- (28) 

i 

in terms of analytical representations of the coefficient functions bi and on 
a domain D = {0,T] x Q with Q C M" a bounded domain. The expansion is 
local in time (as is the WKB-expansion). Our explicit expansion is derived 
from the ansatz 



1 / y;"_ Ax" 

P{t, X, y) = — n exp ' 



-) {^du{t,x,y)t'^^ . (29) 



Here, Axj := (xj — yi). Note again that this is different from the WKB- 
expansion which is of the form 

p(t, X, y) = exp - ^^=\ ' + ^ cfc(i, X, y)t^ . (30) 
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The coefficients dk in ()29p are more difficult to compute than the coefficients 
Cfc in (j30p . However, concerning the growth of the coefficients with respect 
to the spatial variables we expect for coefficients with bounded derivatives 
and fixed time that 

Ck ~ Ax^^ (31) 

holds formally for the WKB expansion and this implies that higher order 
approximations involving for k > 2 may cause problems (we have no 
negative sign for Cfc on the whole domain in general). 

Remark 1.3. In the situations of low volatilities as considered in [2] it is 
sufficient to compute the terms cq and ci even in order to get accurate 
results options with ten years maturity in a scheme with one time step. 
However, in general we need higher order approximations. 

Hence we expect for coefficients with bounded derivatives and for fixed 
time that 

4 ~ exp(co)Ax2'' (32) 

holds formally for expansions of the form ()24p . However, as the growth of 
Co is linear in Ax the highest order term of (|24p is an effective damping 
factor as |Aj;| becomes large (note, however, that we consider only bounded 
domains). Expansions of the form (j24p were considered in a more general 
framework in proofs of the Atiyah-Singer index theorem of course. However, 
in that case only the behavior in the limit t J, is of interest. However, from a 
perspective of computational implementation it is also of interest to consider 
how large the time horizon can be chosen given a certain (bounded) domain 
and a certain set of coefficients bi such that an expansion of the form ()24p 
holds. First we do the analysis in the special case of (|28p . We shall assume 
that for all 1 < i < n the functions (t, x) — t- bi{t, x) are of linear growth and 
equal their Taylor expansion and have bounded derivatives of polynomial 
growth, i.e. 

\D]^D^bi{t,x)\ <C"+^"^ (33) 

for some C > 0. Since we consider bounded domains and from the perspec- 
tive of computation an important case are finite Fourier series, i.e. 

mo 

biit, x)= + + bj cosift + j-x), (34) 

j=-mo 

where uiq = (mQ,mQ,-- - jTUq) and j = )i") are n + 1-tuples, 

and j ■ X := Y17=oji-^i denotes the scalar product of dimension n where we 
identify xq with time t for simplicity of notation. Since all L^-functions can 
be approximated on a bounded domain up to any degree of accuracy the 
class of finite Fourier series is satisfying from a computational or practical 
point of view. In the more general case of variable diffusion coefficients a*. 
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we assume that The result is also fundamental for investigation of the more 
general situation of equations with spatially dependent coefficients 



or equations with coefficients dependent of space and time as in 

|-E4('-)^ + EMM)|^, (30) 

because some arguments will of the reducible case will transfer to the irre- 
ducible case. However, the extension hinges on a deeper analysis of (fT2]) . We 
will state the result but an extension of a the proof is given in a subsequent 
paper. Since second order PDEs are essentially symmetric, it is essentially 
an equation of form 

1 " 

= 4 X] ^i(^)'^R,x,d%^i, (37) 

i=l 

which has to be solved. For an extension of the convergence proof additional 
assumptions on the coefficient functions are needed. In order to compute a 
lower order bound of convergence we assume that 

|L>,"D>y(t,x)| <C"+I°I (38) 

for some C > and for all 1 < i, j < n. Since we consider bounded domains 
and from the perspective of computation an important case are finite Fourier 
series, i.e. 

mo 

aij{t,x)= a* sm{j'^t + j ■ x) + b* cos{ j^t + j ■ x), (39) 

i=-mo 

with an analogous notation as above. In the following we state a main 
theorem which makes the preceding remarks precise and then prove the 
theorem. Then in section 3 we compute a recursively explicit formulas for 
the coefficients d^- In section 4 we prove convergence and determine a lower 
bound for the time horizon where the expansion converges. In section 5 we 
consider extensions to parabolic equations with time-dependent coefficients. 
In section 6 we consider the expansion of parabolic equations with variable 
coefficients. 

In section 7 we consider application to (micro)-hypoelliptic and semi- 
elliptic equations, i.e. the design of weak higher order schemes in this con- 
text. In section 8 we consider applications to American derivatives as a 
second type of semi-elliptic equations which are not (micro)-hypoelliptic in 
general. 
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2 A convergence result for reducible diffusion equa- 
tions 



While heat-kernel expansions of the form ()29p are well-known the following 
questions deserve investigation: 

• what is a lower bound for the exact time horizon Tq where such an 
expansion holds ? Can Tq be computed in terms of the coefficient 
functions 6j? 

• is there an exact formula in terms of analytical expansions of the co- 
efficient functions? 

Both questions are very important as they are fundamental in order to com- 
pute an efficient scheme for parabolic problems of type ()28p . The answer to 
the first questions tells us what time step size we have to choose in order 
to construct a locally analytic weak higher order scheme. The answer to 
the second question provides a formula for each time iteration step of the 
scheme, where the semi-group property is invoked to get a global scheme. 
We have 

Theorem 2.1. Given assumption Ii33\) there exists a finite time horizon Tq 
such that on the domain Q x (0, Tq] for any finite Tq > and any domain 
Q C M*^ a constant (3 can he computed such that the fundamental solution of 

ou o u sr^ , ou , , , 

^ = Ea^+E^^^a;^ (40) 

has the pointwise valid representation 

p{t,x,0,y) = -=l=exp (f2dk{t,x,y)tA , (41) 

for j = 1, • • • ,n, and for {t,x) G (0,/3To) x 0,. If jg^D holds then a lower 
bound of the constant f3 is given by 



^ 3(n(2|mo| + l))ei?2|r72o|2' ^ ' 

where 2|mo| +1 is (an upper bound of) the number of terms in the fi- 
nite Fourier representation of bi (any i € {I,-- - along with |mo| : = 
maXjg|o, - n} "^0 '^'^'^ R is a radius of a ball Br{Q) such that the spatial part 
of the domain is included, i.e. -6^(0) 5 and e is an upper bound for 
the Fourier coefficients of the Fourier representation of the drift function bi 
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where i G {I,-"" j'T-}- For the coefficient functions the following holds: 
for k = we have 



do{t, X, y) = exp (^^{Vm - x,n) J b^it, y + s{x - y))dsj , 
dm{t,x,y) = y —drn-k / Rk-i{t,y + s{x-y),y)s''ds 



with 



Rk-iit,x,y) = f Cfc_i + Acfc-i + E?=i Er=d (^c^^Cfc- 



l-r 



(43) 
(44) 

(45) 



Remark 2.2. Note that the relation p36|) contains also the relation of the 
time horizon of the convergence to the size of the domain, i.e. the time 
horizon is proportional to the inverse of the square of the domain. 



Theorem 2.3. More explicitly, we have 

co(t, X, y) = co(x, y)= -Y,i E7 ^(y) Ax^+i^ 

= E7 C07Ax'^ 

and, given the power series representation 

Ck-i{t,x,y) = ^C(fc„i)^i(y)Ax'^t' 

we have 

Ck{t,x,y) = Y.^^ilc(^k-iYii{y)^x^t^+ 

E7 { Ei Ep+a=7(Pi + + l)Cr(/3+li)C(fc-l-r)(a+li) 

+ Ei(7i + 2)(7i + i)cfc(7+2o + Ep+a=7(E i^Hy)^ 

{ai + l)c(fc_i)(„+i^)} (ELo^'fcl^^'') ' 

n 

where with 5s := ^^(^j; and 

/3 



(46) 



(47) 



(48) 



i=l 



6=0 



Sit 



5=0 



(1 - r)5E + k 



(49) 



n 



,(7-5) 



Ax^ 
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3 Formal computation of solution of the parabolic 
equation 

First we consider the equation ()28p with time-homogeneous coefficients, i.e. 
where the coefficient functions x — t- bi{x) depend only on the spatial variable 
X. It turns out that the ansatz 

pit,x,y) = — L^exp ( - ^^=; — ^ ) ( \ ]dk{x,y)t' ) (50) 



leads to first order partial differetial equations for with variable coeffi- 
cients which are difficult to solve if looked at from an abstract point of view. 
Hence we compute the dk via the Ck of the WKB- ansatz 

Pit, X, y) = -= exp - ^'=1 ' + ckix, y)t^ . (51) 
"^^^^ V fc=0 / 

Hence we first derive recursive relations for the coefficients Ck and then 
get the recursive relations for the dk via the general logarithmic recursion 
outlined in the introduction, i.e. we use the recursion 

(io = exp(co), (52) 

and 

k . 

dk = 'Y^dk-iCi. (53) 

4 = 1 

In a second step, using global analyticity of the coefficient functions bi, we 
derive the explicit solution in terms of Taylor power series of bi. For the 
time derivative we get 

For the first and second spatial derivatives we get 
dp I —Ax I 



dxi \ 2t 



(55) 



and 



' -h + Ek£ck{x,y)t^ 



dxf I 2i ' ^fe dx^ 



(56) 



+ {-^ + j:kI-Mx,y)t^y]p{t,x,y). 
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Plugging into ()28p and ordering with respect to the terms t ^ etc. we 
get the following recursive relations for the cj^, where 1 < j < n: 



and for all /c — 1 > 0. 

+ Em bmix)g§^Ck-l = Rk-l{x, y). 

(59) 

Note that the first order coupling of the system is essentially reflected in the 
recursive first order partial differential equations starting from (|59p . This 
would be different if we had coupling via the second order terms and it makes 
the solution of the system much easier. Note that equation ()57p is satisfied. 
Equation ([58|) is equivalent to 

^ ~ E bm{x)Axm, (60) 

I ^' m 

with the solution 

co{x,y) = S2iym- Xm) 6; (y + s(x - y))ds (61) 

■^0 I 

and for all k > 1 we have 

Ck{x,y)= Rk-i{y + s{x - y),y)s^ds (62) 
Jo 



with Rk-i as in equation (|86p . Next we compute the solution explicitly 
doing the integral for cq first. We abbreviate Ax = {x — y) with components 
Axj = (x — y)i and for a multiindex a = (ai, • • • ,«„) we write Ax° := 
n"^-i^Ax"\ Furthermore, we define |a| = Yli^i 

bm{x) = Y,-^h^,^{y){/\x)\ (63) 

7 ^' 
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along with some multiindex 7, then 

Co{x, y) = - Ylm ^^rn /q^ E« ^rrXv + sb.x)ds 



^"m-f\y )'-^-^ 1+171' 
1 

i+ItI 



Next we compute for A; > 1. We have 

+Acfc_i + Em ^'m^^Ky + s{x - y))s''-^ds. 
Assuming that Ck-i equals its Taylor series for every y £ EJ^, i.e. 

Cfc_i(x) = ^C(fe_i)^(y)Ax'^, 

7 

then we may evaluate the derivatives occurring in Rk^i as follows 
-Q— = 2^(7i + l)c(fc_i)(^+i,)(y)Aa;^, 

* 7 

and 

= 5^(7. + 2)(7. + l)c,(,+2,)(?/)Ax^ 

and 

|^^^|^ = J^| ^ (ft + !)(«, + l)c,(^+i,)C(fc_i_^)(«+i^) 

* * 7 [/3+a=7 
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For the multiindex 7, we have 



{y + s{x - y)y's^-^ds 



(70) 



/o 



HE 

1=1 y<5i=o 



,("i-'5i)^^5i^<5j 1 ^fe-1 



7i 



^ 1 



<5s + /c 

I 



5=0 

7 



7i 



n V<5.!(7.-5i)! 



.,(7-5) 



<5=0 



where 5s := E*^' ^'^'^ ~ ^=\^^^ — • Hence 



1=1 



1^7 { YLi J2l3+a=-yil^i + + l)'>(/3+l .)C(fc_l_r-)(a+li) 

+ Ei(7i + 2)(7i + l)Cfc(7+2,) + E/3+a=7(E wMmAy^^ 



(71) 



(ai + l)c(fc_i)(„+i^)} {EhoPfs^^^) ■ 



4 Proof of convergence and computation of the 
time horizon of convergence 

It is essential to show that for fixed x,y £ il. (where il. is some bounded 
domain) 



k=0 



(72) 



represents some analytic function on some interval [0, /3Tq] for some (3 to be 
determined. Then it follows that 



(73) 



k=0 



Kk=0 



is also an analytic function on that interval. 
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Since C is bounded, there is a ball Bji{0) around with radius R 
such that Q C Br{0). Recall that 

coix, y) = 'y2iyra - Xm) / &m(y + s{x - y))ds, (74) 

and for all > 1 we have 

Ckix,y)= Rk-i{y + s{x -y),y)s''~'^ds (75) 
Jo 

with 

Rk-iit, X, y) = Acfc„i + X:r=i Eto (^c^^cfc 



(76) 



Let us first consider cq. The coefficients of equation (|34p have a representa- 
tion of the form 



mo 

6i(t, x) = ^ exp{i{ft + j ■ x)) + fj exp{-i{ft + j • x)), (77) 

i=-mo 

for some real numbers ej and fj. Let 

e = max {ej,fj\jeJ} (78) 

where J is the set of all multiindices of the sum in ()77p . Then we get the 
estimate 

Coix,y) = YZi=l{ym - Xm) jlT.mbm{y + S{X - y))ds 

(79) 

< ei?(2mo + 1) =: Cq. 
Next the time transformation 

t = pT (80) 

transforms the equation 

j = l J k 

into the equation 

j=l J k 
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where u{t,x) = v{t,x), and where ^ = ff|^ = ff^- The analogous 
representation of the solution is of the form 



1 / V^*^ A 2 



Vi^" V 4/3t 



■) + ' (83) 



for J = 1, • • • , n. Plugging ()83p into ()82p and collecting the terms with r ^, 
etc. we get (we feel free to write t instead of /3t if convenient) 



r 



4/3t^ ' ^4/3V 



and for all A; — 1 > 

r'^-i: kc,,^ + ^Ei^^i'-^= /3Ac,_i,^ + /3Er=iE'=o'(4c.,^ic,_i 

+/3 Em ^m(^)a£:Cfc_i,^ =: Rl_^{x,y). 

(86) 

We divide equation (f86|) by /? and get the solutions (the solution for cq,/? 
equals exactly that for cq in ([71|) ) 

CkA^^y) = -Rfc_i(y + s(a; - y),y)s'3 ds. (87) 

JO 

Next we prove 

Theorem 4.1. For a give time horizon Tq > there exists f3 > such that 

for each x,y £ Q,l < j < n Ck,(s{x, y) i as k '[ oo. (88) 

Proof. Next, a majorant of Ck^pix^y) is obtained as follows: we consider 
three types of operators O^'"", O^' with positive integers fc, and acting 
on a single function / :J7xr2— )-Morona families of functions {fi)i<i<k ■ 
f2 X — )• M, namely 

Ofc'" [/] {x,y) := lAf{x,y) 

d"^ [/., •••,/!] (X, ,) := f Er=i Eto t %f (89) 
Of [/] (^,2/) :=|Em^n.(^)at/(^,2/)- 
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Let 
For 

we have 
where 



Mfc :={(afc,-- - ,«i)|ai G {1,2,3}} (90) 



'^M •= sup Ck,ii{x,y) (91) 

' x,y(iQ, 

cl^f}< sup ^ O°'"co(x,y). (92) 



O^'- [/] (X, y) := Of '"O,"--'^ o . . . o or-" [/] (rr, y). (93) 

First let Ifc (resp. 2^,3^) the multiindex a S such that for each 
1 < m < k am = 1 (resp. Om = 2, = 3). Hence 

Ol [f] ix,y) = (Oi)' [/] ix,y) = A^' [/] {x,y) (94) 

etc.. In order to compute this expression we first apply the Laplacian to 
exp(moix). We get 

|A''exp(imoa;)| < n^lmol^'', (95) 

where |mo| := niaXjg|o, i, ■■■.«} '^o- Since there are n drift functions bi the 
explicit representation of cq has n(2|mo| + 1) terms with the factor exp{i{j ■ 
x)). Since the dimension is n we have n factors 6xi := (xj — yi) in the 
definition of cq. Note that for m > 



^) (feiexp(i(j -x))) 



(96) 



= i&jj exp(i(j • x)) + 5xi l^-^j exp(i(j • x)). 

Hence, 

- sup |A%(x,y)| <n(2|mo| + l)e(fc + i?) ^ ^ ' °' ^ 10 (97) 
as /c t if 

< /3 < (98) 
n|mo| 

Next we observe that on a domain of radius 

l^3,nr , ,w ,|^-fc^ /?'fc(fc!n^(2ko| + l)fe) , ^ 

[co(x,y)] (x,y)| < e ^^0 (99) 

as /c t CO if (note the strict inequality sign because of additional factor k) 

< /3 < , \ — . (100) 

e(n(2|mo| + 1)) ^ ' 
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The operators of quadratic type applied to co{x,y) 0^'"co decrease also to 
zero as /c t c« if /3 is small. Inductively with respect to k you show that 

|0f [exp(imox)] (x, y)\ < |moP'=3'=A;!. (101) 

We get 

l|^2,nr w „ ^ /3^(n(2|mo| + l))^#i?2^|mop^3^^fc! , ^ 

[coJ(x,y)|< 10 (102) 

for /c t oo if 

^ ^ 3(n(2|mo| + l))ei?2|mo|2 ^^^^^ 

For large A; this is essentially the largest term of all the 3^^ contributions 
in the sum (j92p for large k [k fixed). You can show that if (jl03p holds, then 
we have for k > ko (some ko > 

|3^0gM(x..)l< ^^^^'=^^''"'f^'='"°'"' iO. (104, 

as fc t oo, and this is also the estimate which holds for Ck for large k. Here 
we choose /3 such that in a summand in O^'"co(a;, y) in (|92|) each occurrence 
of an operator of form O^'" can be replaced by an operator of form O^'" 



in order to get a majorant estimation. So in the sum in (j92p it suffices 
to concentrate on the summands consisting of concatenations of operators 
of form O^.'" and O^.'"- For natural numbers / let us define an increasing 
sequence of numbers ki < k2 < • • • < ki < ki^i • • • , and operators 



01,n /ol-" /^Ij" 

, , := U, o ■ ■ ■ o (J, 
ki+ih fci+1 ki 



02,n /^2.n /^2,n 

, , := U, o ■ ■ ■ o (J, 



(105) 



Then in the summands o (|92p we have to consider the asymptotic behavior 
of values of family of operators of form 

Oj^.fc, o Ofcife,^, ° • • • ° ° (106) 

or of form 

02,n ,^1," /-ill" /o2,n /in^yN 

fci+ifc, ° o . . . o O,;,^ o O^;;,^ (107) 

applied to co(a;, y) as A; f oo. If there is only a finite occurrence of operators 
of form O^'" in such a family f ()106p or ()107p ). then the asymptotic behavior 
is clearly the same as for O|co(x,y). If on the other hand there are infinite 
occurrences of operators of form 0]j^ in f ()106p of ()107p ). then for large k 
Of,co{x,y) becomes a majorant of such a summand. Hence, the estimate 
p04p is a majorant for large k and proves the convergence of the series in 
(1921). □ 
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5 Generalisation to the time-inhomogeneous case 



We start with the formal computation in the time-inhomogeneous case and 
then show how the local convergence proof of the previous section can be 
extended. 



5.1 Formal computation of recursive coefficients in the time- 
inhomogeneous case 

We consider parabolic equations with time-dependent coefficients of the form 
- + A. + ^5.(t,x)— = (108) 

We consider the ansatz 

1 / Ax^ ^ \ 

p{t, X, 0, y) = ^= exp — + V Cfc(t, x, y)t^ . (109) 

Compared to the time-homogeneous case the time derivative contains an 
additional term. We have 

^{t,x,0,y)=(^-% + ^ + Er=o ^ (t, X, y)t^ 

+ Hk ^cfc(t, X, y)t''"^ j p{t, X, 0, y) 

(110) 

The spatial derivatives are essentially the same as in the time-homogeneous 
case. We compute 



+ X] -^^kit, X, y)t'' j p{t, X, 0, y), 



^(i, x,y)=l + > J £rCfc(t, X, y)e ) p{t, x, 0, y), (111) 



and 



(112) 

+ (-^ + Efc ^Cfc(t, X, y)tk) ' j Pit, X, 0, y). 

Plugging into equation (|28|) and ordering with respect to the terms t~'^,t~^ 
etc. we get the following recursive relations for the c^, where I < j < n. 
First, the highest order terms are the same as before: 

I 
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The terms of order t ^ are essentially as before (we just have to add the 
i-argument in the coefficient functions bj): 



i = -E^-ifeA^-'g-EMM)Aj. (114) 

/ \ I ' m / 



For /c — 1 > we get an additional i-derivative on the right side: 



+ Em ^rn{t, x)^Cfc_i = i?fc_l(x, y) 

Hence, 

\ " At-, 

dxi 



(115) 

^ Ax;^ = -y2bm{t,x)AXm, (116) 



I 

which has the solution 



co{x,y) = '^iym- Xm) / '^bi{t,y + s{x - y))ds, 

m I 



(117) 



and for all /c > 1 we have 

'1 



Ckix,y)= Rk-iit,y + s{x-y),y)s^ds (118) 

JO 

with Rk-i as in equation (jllSp . The explicit calculation of the solution is 
know completely analogous, so it suffices to write down the results. We write 

bm{t, x) = Y^ ^b^,^it, y)iAxy (119) 

7 ^' 



(120) 



along with some multiindex 7. Then 

co(t, x,y) =- 6,„^(y) Ax^+i' 

= E7C07(t,y)Ax'^ 
Given that Ck-i equals its Taylor series for every y G M", i.e. 

Ck-i{t,x) = ^C(fc_i)^(t,y)Ax'^ = ^C(fc„i)^;(y)Ax^t^ (121) 

7 7:' 
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we have 



Ck{t,x,y) = E7,/^C(fc_i)^/(y)Ax^t' 



where the are defined exactly as before. The proof of convergence is 
then analogue to the time-inhomogeneous case. 

6 Generalisation in the case of variable coefficients 

In order to provide a first discussion to various types of semi-elliptic equa- 
tions let us formally compute the local expansion 



p{5t,x,y) = ^=^exp + Y,Ck{x,y)t' , (123) 



(122) 




of the parabolic equation 



du 
'dt 



ij i 




(124) 



First we compute for t > 



dp 

dt ~ 



( 



i + + EZoik + l)cfc+i(x, y)t'') p, 





2t ^ l^k=0 dxi V-^' y)'' 



) 




(125) 
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Plugging into (|124p we get 

( - i + + EZoik + l)c.+i(x,y)t'= - i E,, <,(^)x 



^00 dc, 



2t + I^fc=o ax. 



oo 9c, 

fc=o ax 



— (126) 



( - i + %^ + Er=o(fc + l)cfcH-i(x, 



2i 



oo acfc^fc 



"I" Z^fc=o 1 2^i=o axi dx 



2t 



(127) 



&.(x) ^- 1 + Er=o y)t') )p = 0. 

Collecting terms of order we have 



(128) 



Note that here d^. is the partial derivative of cP with respect to Xj. Equa- 
tion ()128p is closely connected to a Hamilton-Jacobi equation and admits a 
unique solution if the boundary condition, i.e. the condition d{x, y) = if 
X = y, is satisfied. Collecting terms of order we get 



- ^ + hd' + - y a*,(x)f^|^(x,y) + ^^{x,y)) = 0. (129) 
Equation (|129|) is a linear first order equation which can be written as 

,y) = o. (130) 



n 1 
2 



As we shall see later, this equation together with some boundary condition 
co(y,y) = -^ln ^det(a«i(y)) (131) 
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determines cq uniquely for each y € M". In general, the boundary condition 
on Co ensures that p is a probability density (i.e. integrates to 1). This is not 
essential as far as existence, uniqueness, and convergence of the coefficient 
functions Cfc are concerned. If we define co(x, y) — co{y,y) =: co{x,y) then 
Co satisfies the equation ()130p too with the boundary condition co{x,y) = 
if X = y. For A; + 1 > 1 we get 



(k + l)c,+i(x, y) + i E., <,ix) + 

(132) 

with the boundary conditions 

Ck+i{x, y) = Rk{y, y) if x = y, (133) 

Rk being the right side of ()132p . 
We have 

Theorem 6.1. // the assumptions (3E\) o-nd / fggj) are satisfied, then there 
exists a finite time horizon Tq such that on the domain x (0, To] for any 
finite To > and any domain Q C M"- a constant j3 can he computed such 
that the fundamental solution of 

du d'^u , du 

has the pointwise valid representation 



1 

p[t, X, 0, y) = n exp 

\/4vrt(r) 



■^%^) (T.^dk{t,x,y)t'^^, (135) 



for j = 1,--- , n, and /or (t, x) G (0, /JTq) x ^2. // and (E^j /ioZd i/ien a 
lower bound of the constant /3 is given by 

^ ^ 6(n(2|mo| + l))efl2|rno|2' ^^^^^ 

where 2|mo| +1 is (an upper bound of) the number of terms in the fi- 
nite Fourier representation of bi (any i £ {1, • • • ,n} along with |mo| : = 
niaXjg|o^...n} 77T-Q and R is a radius of a ball -B_r(0) such that the spatial part 
of the domain 0, is included, i.e. Bji{0) 5 ande is now an upper bound 
for the Fourier coefficients of the Fourier representation of the drift function 
bi where i S {!,••• ,n} and of the Fourier coefficients of the Fourier rep- 
resentation of the diffusion functions a-ij, 1 < i, j < n. For the coefficient 
functions dk the following holds: for k = we have 

do{x,y) = exp{co{x,y)) , (137) 
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k 



dmit, x,y) ='y —dm-kCk, (138) 
where Ck is defined as above. 



m 

k=l 



7 Application to linear semi-elliptic equations 

First we consider the application to (micro)-hypoelliptic equations. Then 
we consider the appUcation to semi-ehiptic equations which are (micro)- 
hypoehiptic on a hnear subspace. More details about the analysis of semi- 
linear equations will be provided in version 3 of (which will appear in 
arXiv shortly). 

7.1 Application to (micro)-hypoelliptic equations 

In the situation of linear semi-elliptic equations described in the intro- 
duction consider the case n = d, i.e. consider a matrix-valued function 
X — )■ {aji)'^'"^{x), < i,j < n on M", and m + 1 smooth vector fields of 
dimension d 

d 

^« = ^%i^> l<i<m, (139) 
j=i 

where < i < ?n. Consider the Cauchy problem on [0, T\ x with time 
horizon T > 0: 

t = h Er=i A'" + Aon 

(140) 

u{Q,x) = fix). 

Assume that (|14U|) satisfies the Hormander condition with respect to the 
subspace W^, i.e. assume that 

{A,[Aj,Ak],[[Aj,Ak],Ai],---\l<i<m, < j,k,l ■ ■ ■ < m} (141) 

spans at each point x. The existence of regular solutions of the Cauchy 
problem in this case is well known (cf. [8] ). Indeed Hormander's result 
shows us that there exists a family of smooth transition densities if (jl4ip 
holds for every x G M'^ . Indeed we may extract the following Malliavin 
estimate of the density from |17| . 

Theorem 7.1. Consider an n- dimensional diffusion process associated to 
{l-jC^ of the form 

n d 

dXt = HXt)dt + (Jij{Xt)dWi (142) 
i=i j=i 
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with X{0) = X £ M" with values in M" and on a time interval [0, T], i.e. 
assume that the solution of the Cauchy problem has the probabilistic 

representation 

u{t,x) = E^f{Xt)) (143) 

Assume that bi,aij G Ci^- Then the law of the process X is absolutely 
continuous with respect to the Lebesgue measure and the density p exists and 
is smooth, i.e. 

p : (0, T] X R" X M" ^ M G C~ ((0, T] x M" x M") . (144) 

Moreover, for each nonnegative natural number j, and multiindices q, (3 there 
are increasing functions of time 

^j,a,/3,5j-,,;3: [0,r] ^M, (145) 

and functions 

nj,a,p, mj^a,f5 : N X N" X N" ^ N, (146) 

such that 

\DiD^D^^pit,x,y)\ < ^^^^^^4S^^^P (147) 

Moreover, all functions |j^5[ ) and depend on the level of iteration of 

Lie-bracket iteration at which the Hormander condition becomes true. 

Now consider the recursion equations (|129p and ()132p for cq and Ck+i for 
A: > of the expansion (|123p of the last section. 

Let / : — )• R be a continuous function of atmost exponential growth, 

i.e. 

for all X G M" |/(a;)| < Cexp(C|x|) for some constantC > 0. 

Assume that Ojj , 6j are bounded coefficients with bounded derivatives with 

|-D>ij| < \Dl^aij\ < Cl^l (148) 

for some C (e.g. aij and bi have representations in form of finite Fourier 
series). Then under these conditions using the [?] one may prove that 

lim/ f{y)pDf^(t,x,y)dy= f{y)p{t,x,y)dy, (149) 

where p denotes the density of the process X according to theorem (|7.ip 
above. Then we may approximate the density via expansion on D^^^ where 
it is strictly elliptic and set up an higher order scheme. Details of analysis 
will be given in a subsequent paper. Note that this observation can be used 
to construct weak higher order estimates for diffusion market models which 
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satisfy the Hormander condition. Recall that an approximation scheme Y 
converges weakly with order 7 > to X as At 4, and with respect to a 
function class C, if for all g G C 

\E igiXr)) - E {giYr)) \ < CAT"' (150) 

as AT 4, 0. It depends on the regularity of functions in C whether ()150p is 
a strong condition. In finance we may low regularity of payoffs, e.g. in the 
case of digital payoffs. This leads to sophisticated weighted Monte-Carlo 
schemes of bounded variance. The related estimators in [2] may be adapted 
to the more general situation in a quite straightforward way (well, the proof 
of bounded variance and the error estimates become a little more intricate) . 

7.2 Applications to semi-linear equations which are hypoel- 
liptic on some linear subspace 

Note that for n > d the Cauchy problem (|140p does not satisfy the Hor- 
mander condition on the whole space M" in general but only on a linear 
subspace. Especially, a density may not exist in a regular sense. As an 
example, consider the following Cauchy problem on [0,T] x M?: 

(151) 

^(0, X) = f{xi, ■■■ ,Xd)+ g{xd+l, • • • , Xn). 

Assume that the Hormander condition is satisfied on a subspace of dimension 
d (for any (x^+i, • • • ,Xn) fixed). Hence the fundamental solution of 

du 1 7., , , d'^u ■J^ du 

i,j=l ■' j=a+l 

exists (for any (X(i+i, • • • , x„) fixed), and the solution of (jl5ip becomes The 
solution of this equation is 

/(yi, • • • , yd)Pd{t, X, y)dy + g{xd+i + Hd+it, ■■■ ,Xn + fJ^nt). (153) 
This leads us to a 'distributional density' of the form 

p{t, X, y) := pd{t, X, y)6{xd+i + /^d+i* - Vd+i, ■■■ ,Xn + l^nt - yn)- (154) 

We see from this example that a density exists only in a distributional sense. 
However, if the initial data / satisfy an exponential growth condition, and 
are smooth on the space and locally L'p on the subspace M'^ where the 

Hormander condition holds, then the Cauchy problem ()140p has a regular 
solution, i.e. a solution in C°° ((0, T] x M"). The simple example shows that 
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in general there exists no regular density in a situation of degenerate diffusion 
models with n > d. Especially, the regularity theory of densities in the 
context of Malliavin calculus does not apply directly (as we remarked in [7] ) . 
However, the regularity theory for densities of Malliavin calculus may still 
be useful, especially for the analysis on the subspace of dimension d where 
the Hormander condition holds. Let us consider ()140p from a probabilistic 
perspective. The associated Stratonovic integral of a process starting at 
x&W is 



where Wi denote Brownian motions and o indicates that the integral is in 
the Stratonovic sense, li n = d and the Hormander condition holds, then 
the associated covariance matrix process at is almost surely invertible, i.e. 



for every real number p and t in some arbitrary finite time horizon [0, T]. 
Clearly this is not true in the example above. Note, however, that the asso- 
ciated process exists, i.e. global existence for ordinary stochastic differential 
equations as is well-known in the context of elementary stochastic analysis. 
A standard theorem of ordinary stochastic differential equations (for state- 
ment and proof cf. [12]) shows the existence of Levy-continuous solutions. 
Next we state and an extension of a theorem in [7] (the proof will be given 
in the third version of [7] ). One can use estimates obtained in [T7] which 
we cited above. In order to see how analytic expansions can be used in 
this context, our main interest here are some of the constructive aspects of 
the scheme which leads to the global existence and regularity proof. Let's 
consider the theorem first. Our interst here is the use of analytic expansions 
in this context. A detailed proof will be given in version 3 of [7] in arXiv 
shortly. Let us consider the time-homogeneous case. Consider a matrix 
function x — )• {vji)'^'"^{x), 1 < i < ra, < z < m on M^, and m smooth 
vector fields 



where < i < m. Consider the Cauchy problem on [0, T] x M" (where T > 
is an arbitrary finite horizon) 



Xt = x + 




(156) 




(157) 



/ du 



1 V'" T/2 



U + VqU 



dt ~ 



< 



(158) 
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This may be rewritten in the form 



(159) 

u(0,x) = /(x), 

where 

m 

fc=i 

A general reduced Cauchy problem is defined by the condition that for all 
t G [0, T] and x G M" 

«,) (^) (161) 

has rank d = d{x) < n, where for each x £ M"" the number d{x) is determined 
by the Hormander condition at x € M", i.e. 

{y,(x), [Vj, Vk] {t, x), [[Vj, Vk],Vi]{t,x),---\l<i<m, 0<j,k,l--- <m} 

(162) 

spans a linear subspace Wx of dimension d{x). In this case we may consider 
the intersection of all x-dependent subspaces which are spanned by the local 
Hormander condition at x, i.e. 

Ih ■■= n^m-Wx (163) 

where Wx is the vector subspace of dimension d{x) spanned by ()162p above 
at X. Here the notation Ih indicates that we consider a intersection of spaces 
defined by local Hormander conditions. Our most general theorem will show 
that regular global solutions of ()158p (rsp. ()159p ) exist if the data are rough 
(i.e. in (M") only in Ih and are smooth on the complementary vector 
subspace M" \ Ih- In the following section we look at the situation from 
the perspective of elementary stochastic analysis. We observe that from this 
perspective regularity is closely linked to regularity of the data. We shall 
see then that a constructive analytic scheme together with Malliavin type 
estimates lead us to stronger results. Next consider a matrix- valued function 
X — )• (uji)"'"*(x), 1 < j < n, < i < m on M", and m smooth vector fields 

V, = 5^t;,.(x) — , (164) 

where < i < m. 
We have 
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Theorem 7.2. Let I < p < oo. Consider the Cauchy problem Iil58\) on 
[0,T] X M". Assume that the initial data function f : — )• M satisfies 



(i) i/ie function f is Lf^^, 1 < p < oo on Ih, 

(a) the function f is C°° on W^\Ih, 

(Hi) for all x G 

|/(x)| < C exp(C|x|) for some constant C > 0. 



(165) 



Assume that the coefficients are smooth (i.e. C°°) of linear growth with 
bounded derivatives, i.e. 

Vji e Cfl (M") (166) 

for i = and 1 < j < n, or I < i < m and 1 < j < n. Then the Cauchy 
problem II158\) has a global classical solution u, where 

u G C°°((0,T] X M"), (167) 

where the singular behaviour is t = is determined by the Malliavin-type 
estimate in ]j7| / as follows: for given natural numbers m and N there is a 
number q such that the solution u and its time derivatives up to order m 
and its spatial derivatives up to order N are located in the space. 

Cl^N ([0' T] X M") := {v\t'iv G C^,n ([0, T] x M")} , (168) 

where 

C^,^([0,r]xM"):=i/| + mf\\<^ 

l<m \ol\<N 

(169) 

and \\.\\ denoting the supremum norm. Moreover, q = max^a\<N '>T'm,a,o — n/2 
where n^.^afi is determined by the estimate in fF^ of the singular behavior 
of the density. 

Since we are interested in the apphcation of analytic expansions let us 
consider the main steps of the proof which lead to the constructive scheme 
(cf. [7] for details). For most applications in finance the situation in ()7.2p 
can be reduced to the block structure, where the whole space M" can be 
decomposed into a part where the Hormander condition holds and a 
complementary part W^~'^ where the operator looks like a vector field. For 
associated first order equations with source term one observes (for a proof 
cf. 0) : 

Proposition 7.3. Fix G M"'. Assume that the conditions of theorem 1 
are satisfied. Assume that g £ ([0)^] ^ W^). Then there exists a smooth 
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global flow 3"* generated by the vector field below on [0, T] x M" such that 
the first order equation problem 



(170) 



has the solution 

u{t, x"^, x"-'^) = / (^x"^, + g{s, x"^, 3^-'x''-'^)ds. (171) 

The notation above which indicates that some coordinates are fixed {x'^ 
or x"'~'^) is a little cumbersome and we shall drop it sometimes writing 
just X instead of (x'^,x"~'^) in the following when it is quite clear from the 
context which components of x should be considered to be fixed. Next we 
define a local iteration scheme involving global flows of the type discussed 
in Proposition 17.31 above and solutions of parabolic equations of form 

i,j = l •' i=l 

with x""*^ fixed. The natural ansatz is an AD-scheme of the following form: 
we define 

Vector Field Step: [l > 0) 



at ^i=d+i h^t\-^ J dx. 



u=d^ii^d^y- 



Et=i4(-)S^ + Eti^.(-)^ if ^/O (173) 



i,j=l"'ij\-^J dx.dxj ^ /^i=l J dx 

if / = 



and 

DifTusion Step: (/ > 1 



dt ^i,j=l°'ij(^) dxidxj I]j=lMj(^) dxi 



(174) 



For each m we define u'"(0, .) = /(.) and u"^~^^{0, .) = /(.). Here, in equation 
()174p we understand (x^+i, • • • , x„) to be fixed, and in ()173p we understand 
(xi, • • • ,Xrf) to be fixed. In order to prove convergence time step by time 
step we rewrite the scheme in time-dilatation coordinates {p will be small) 

t : [0,oo) [0,oo), 

(175) 

t(T) = pT. 
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Then we get an equivalent equation in r where the coefficients of the symbol 
of the operator become small if p is small. We have 



dT 

and 



dt 



(176) 

n(0,x) = /(x). 

An iteration step of the scheme considered in transformed time r for 
some time horizon [0, Tq] is then given by 



dr Z^i=d+1 fMiI-^J dxi 

and 

St P l^i,j = l"'ij\-^> dxidxj l^i=l PPn-^J dxi 



(177) 



(178) 



2!-2 



for Z > 1. We start the scheme with 

^-Er=.+im(x)%^ = 0. (179) 
The initial conditions are 

U'''"'{0,X) = f{x), 771 >0, (180) 

where for m = 1, 3, • • • {xd+i, • • • , Xn) is fixed, and for m = 0, 2, • • • (xi, • • • , x^) 
is fixed. The solution can be constructed in the form in the form 

u^ir, x) = uP'\t, x)+Y1 SuP''^^+\t, x), (181) 
i>i 

where for I > 1 

SuP'^'+^ = w^'^'+i - u^'^'-i (182) 

satisfies 

dr P 2^i,j=l"'iJ\-^) dxidXj 2^i=lPPiV'\J>^-^) dxi 

= T.l=d+,PPr{t{rU)^-^, 

(183) 
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and in each substep where the right side in ()183p 
satisfies 

(185) 

Moreover, for m > 1 6uP'"^ has zero initial conditions, i.e. 5u^'^'^''{0,x) = 0. 
For small p the scheme just described is locally convergent with respect to 
time. Then iteration of the scheme in time using the semigroup property 
leads to a convergent scheme of a global solution to the Cauchy problem 

(186) 

uP{0,x) = f{x). 

Note that iteration in time means that we start the next time step with 
the initial data uP(Tq, .), and after repeating the scheme above we get the 
next initial data uP{2Tq, .) and so on. The choice of p depends on certain 
a priori estimates in[T7]. At each time step approximations of the densities 
of the diffusion substeps can be constructed according to the preceeding 
section. This leads to efficient schemes for a considerable class of semi- 
elliptic equations. For problems in higher dimensions probabilistic weighted 
Monte-Carlo schemes are constructed from the scheme above in a natural 
way (we shall discuss them in version 3 of [7] shortly). 



8 Application to American derivatives 

We reconsider the front-fixing method in the case of a multivariate put op- 
tion. We extend the considerations of [3] to (micro)-hypoelliptic operators. 
Our main interest here is how analytic expansions may be used in this con- 
text. Further ananlysis as well as details of global existence and regularity 
proofs can be found in a subsequent paper. We start with the operator 

ij ■' \ i / 

where Vij = {(Ja'^)ij and r may depend on time t and spatial variables S. 
We do not assume that Vij the volatility matrix is strictly elliptic, but we 
assume that the operator is (micro)-hypoelliptic in the continuation region, 
or that the Hormander condition holds in the continuation region. This 
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may be rephrased nicely in the context of the frontfixing method on a half- 
space. Let £ C [0, T] X denote the exercise region and for each t & [0, T] 
let £t denote the t-section of the exercise region, i.e. £j := {x\{t,S) E £}. 
In general beside basic standard assumptions on diffusion market models 
introduced in the previous Section we shall assume that 

(GG) for each t € [0,T] we assume that G £j and that £t is star-shaped 
with respect to 0, i.e. for all S £ £.t and all A G [0, 1] we assume that 
XS e £t. 

Remark 8.1. Note that this means that for a fixed "angle" at S" = 
{Si,-- - ,Sn), i.e. at 



4>s 



So Sr, 



1=1 '-'i l^i=l 



we have one intersection point of the free boundary of the section E-t 
and the ray through which is determined by the angle (ps- 

The condition (GG) is called the global graph condition. The condition 
(GG) holds if X — >• u{t, x) is convex, where {t, x) — )■ u{t, x) denotes the value 
function of an American Put. This is a sufficient (not necessary) condition 
for (GG) to hold. Especially, this condition is satisfied for the multivariate 
Black-Scholes model (Consider the Snell envelope definition in order to verify 
convexity). However, note that from the Snell envelope representation of the 
price of an American index Put 

PA{t,S;K)= sup EqIk-Y,sA (188) 

reStop[(,j,] y j = i J 

we see that the price function pA of an American put option with strike K 
and with maturity T > 

{t,S;K)^PA{t,S;K) (189) 

is homogenous of degree 1 with respect to {S, K), i.e. for any A > we have 
for all t G [0, T] 

PA{t,XS;XK) = XpA{t,S;K). (190) 

Now we can prove 

Proposition 8.2. Assume the American put is written in a market where 
the no-arbitrage condition holds. The price function pA of an American put 
option with strike K and with maturity T > has the property that the 
function 

S^PA{t,S;K) (191) 

is convex for all t £ [0, T] and K fixed. Hence, (GG) is satisfied in this case 
independent of the underlying model. 



34 



Proof. Comparing a portfolio n^(i) consisting of fj, American puts with 
strike Ki and (1 — //) American Puts with strike K2 with a portfolio Pb 
consisting of one American Put with strike K = ^Ki + (1 — /u)/^2 we get 

PA{t,S;K)<f,pA{t,S;K^) + {l-fi)pA{t,S;K2). (192) 

Then for fixed K we may write convexity of pA in S = (/xAi + (1 — l^)^2)K = 
fiSi + (1 — /u)S'2 in the form 



K 



■' (/iAi + (l-MA2) 



(193) 



< ^l\lPA [t, i^; f ) + (1 - ^l)X2PA (t, K- g 



Division of (|193|) by [fiXi + (1 — ^\2) and the observation that pA on the 
left side of ()193p may be rewritten in the form 

/ K \ f K K\ 

p^it.K--——- — =p^(t,K;p— + {l-p)—] (194) 



(/xAi + (1 - ^A2)y V ' ' ^1 ^^2 

along with 

p= , , (195) 

leads to the reduction of convexity with respect to the asset to convexity 
with respect to the strike. □ 

Furthermore, note homogeneity of order one in (5, is a quite natural 
condition. If a selfhnancing portfolio with initial value reduplicates the 
payoff {K — Sr*), then the portfolio Allf reduplicates the payoff {K — St) 

Hence, the free boundary can be written in terms of the angles in form 

{t,qys)^ F{t,cj)s). (196) 
We consider the transformation 

: (0, T) xRl^ (0, T) x [1, 00) x (0, l)''"^ , 

(197) 

'^{ij 'S'l) ■ ■ ■ ) Sn) = ^'j? ) Si ' ' ' ' ' J2'^-i Si ) ■ 

Note that the spatial part of ip ((0, T) x M") is homeomorph to the half 
space H>i = {x £ M"|a;i > 1}. In the following the domain D is the interior 
of the image of V, i-e. D := {0,T) x (l,oo) x (0, 1)""\ We have 



(198) 
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We get 



Ut — Xi g^_^ + 2 Ylij "-ij dx.dxj + I^j dxj + ^ (^1 dxi 

(BCl) u{0, oo, 3:2, •• • , Xn) = on xi = 00 

{BC2) Uxi{t, 1,X2,- ■■ ,Xn) - u{t, 1,X2,- ■■ ,Xn) = -K (-jgg) 

on xi = 1 

{BC3) F{t, X2,--- ,Xn) = K - u{t, 1, X2, • • • , Xn) 

I {IC) u(0, x) = max{K - xi, 0} 

Remark 8.3. We include (BCl) as an implicit boundary condition in order to 
indicate that ()199p is equivalent to an initial-boundary value problem of the 
second type on a finite domain (just by suitable additional transformation 
with respect to the variable xi). 

The mixed condition {BC2) follows from the smooth fit condition to- 
gether with (SC3) . 

Remark 8.4. Note that in the context of market models based on Levy 
processes or, more generally, Feller processes the smooth fit condition does 
not hold in general and one has to be careful concerning generalization at 
this point. 

In order to determine the coefficients af^ and hf we compute first 



TpOXj_ _ ojj -xj uxj _ 1 _ ■sp ( _ \fj_ _d_ _ "^j a (nr\(\\ 
^ dSi ~ xi ' dSi ~ ^ 2^j>2\"i-J Xj) p , QSi ~ dSi dxj ■ \'^^^) 

We observe that 

i i 

It follows that 

EC d dxj Q c dxi d . / o dxi \ d 

i dSi ~ dSi dSi ~ l-^i dSi dxi ^ l^j>2 \2^i dSi J dxj 

= Y^i ^imiwi = Y^i^i{T~ Ylj>2 i^ij ~ ^j) F^j air (202) 



Hence, we have 

^ dSi j ' ""^9x1 
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It is clear that 

ag = 5:.H5,5,^ — , = Y.^kiS,S,^^. (204) 

kl kl 

In order to determine the latter coefficient functions we compute 
Next, for j > 2 we have 

_ % (:|-j-,Sij)({uF + 3:iFi(l-'>u)) 

dxi xiF {xiPy 

Finally, 



(207) 



V i>2 



dSidSk ^ dxi 



dxi 



(209) 



5,iF, + {5,j )F,i {l-Sn)-2FjF, [l-Su ) [5,^ ~Xj ) 



-p^U - Oil) - 2^j>2 W 



Here 5ij is always the Kronecker Delta, Fj is short for Fji is short for 

Q^.Q^^ ■ Now we have determined the explicit form of ()199p . The next step is 
to construct a representation of the solution of ()199p in terms of convolutions 
with the transition density, i.e the fundamental solution related to (|199p . 

Note that only the inner regularity of the free boundary function can be 
proved. Starting with some vP (solution of p99p for F = 1 for example) 
for numerical reasons we may consider an iteration scheme = tvP" with 

= tu^^ and consider for n > 1 the following iteration for 11991 Let 



n-l 

^. Of" I 1 „F"~^ a'v" I ht'"-'^ 

dxj 

(210) 



+r [ -v"-] + u 



where has zero initial condition. On the hyperplane {xi = 1} we require 

<j(t, 1,X2,--- ,Xn) -v''{t,l,X2,--- ,Xn) = -tK. (211) 
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Given u" ^ and F" ^ for each n we look first at the solution for f " in the 
form (note that below contains an integral f^) 



f"(t,x) = /q fj^pyn{t,l,x;T,l,yi)(l)yn{T,l,yi)dHydT 
+ Jd^ u"'~''^{s, y)pi,n {t,x;s, y)dyds, 
where py« is the fundamental solution of 



(212) 



and where solves an integral equation 

i(^^„(t, i,xi) = r^„(t, i,xi)+ 

//o ^H^miP^nit^ ^,xi,T,yi) - py^{t,l,xi;T,l,yi))(t)y„{T, l,yi)dHydT, 

(214) 

(defining (p^^) with 



r»;„(i,a;) = (^-£^py^{t,x;to,y) - py^{t,x;to,yfj t'tljo{y)dy + tK. 

(215) 

Remark 8.5. Asymptotic analysis shows that even VtFt is bounded as t J, 0. 

For the next step we get and via = K—v"'{l,xi) and tu*^ = f". 
In the case of strictly elliptic operators L we get convergence in adapted 
Banach spaces for v"^ (cf. [3]). Regularity resuls and convergence of the 
scheme can be generalised to the case of micro-hypoelliptic operators L, i.e. 
if L has the representation 

^^^EIliA'^ + ^on (216) 

with a matrix-valued function x — t- {aji)'^'"^{x), < i,j < n on M", and 
m + 1 smooth vector fields of dimension d 

d Q 

^i = ^aji — , l<j< m, (217) 
j=i J 

where < i < m such that the Hormander condition with respect to the 
subspace M'^, i.e. assume that 

{Ai, [Aj.Ak] , [[A,,Ak] , A] , • • • |1 < i < "1, < i, /c, / • • • < m} (218) 

spans at each point x G {x\x G M" and xi > 1}. The proof uses the 
estimate in [10] cited above. Details will be given in [3]. 
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